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Abstract 

In the present paper we introduce a new concept of Hardy type space 
naturally defined on the Klein-Dirac quadric. We study different prop- 
erties of the functions belonging to these spaces, in particular boundary 
value problems. We apply these new spaces to polyharmonic interpolation 
and to interpolatory cubature formulas. 

1 Introduction 

In one-dimensional mathematical analysis, Interpolation Theory and Quadra- 
ture formulas are intimately related, cf. |35) , . This relation causes a similar- 
ity between the approaches for estimation of the remainders of Interpolation and 
Quadrature. One approach for estimation of the error in Interpolation theory 
is related to Lagrange formula and uses higher dierivatives of the interpolated 
function (cf. [35], chapter 3, Theorem 4, and [TT], Theorem 3.1.1). The second 
approach uses analyticity of the interpolated function and Hermite formula (cf. 
[35) . chapter 3, Theorem 5, and [11| . Theorema 3.6.1). In a similar way, already 
A. Markov has estimated the error of a quadrature formula for differentiable 
functions in (/) defined on the interval / by means of its A^— th derivative 
(cf. [3S], chapter 7.1, and Davis [TT], p. 344). The second approach estimates 
the error of a quadrature formula for certain classes of functions / which are 
analytic on some open set D in C containing the interval /, (cf. |12j . chapter 
4.6, see also [35], chapter 12.2). However, in both Interpolation and Quadrature, 
the first approach is usually not very practical beyond derivatives of order five, 
see [5^ . 

Interpolation and Approximation of integrals in the multivariate case is a 
much more difficult task. In Numerical Analysis, instead of quadrature formula 
the notion of cubature formula is often used, see [51], [55], [52] and the recent 
survey [lOj . In contrast to the univariate case there is no satisfactory error 
analysis available in the multivariate case, cf. [55| . [S], and part 4 in the last 
Russian edition of the classical monograph [35] . Let us mention that the area 
of quadrature domains which has received a lot of interest recently presents an 
interesting multidimensional alternative and we refer to |14) . 

The present research continues the study of estimates of polyharmonic inter- 
polation and polyharmonic interpolatory cubature formulas initiated in [30j : in 
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the present paper we consider interpolation and cubature formulas in the ball 
in M.'^, while in [3D] the case of an annular region was considered. 

1.1 Gauss- Almansi formula 

In |18| polyharmonic interpolation has been considered for functions defined in 
the ball in Mf^. In the same spirit, in [37] and [2S] we have introduced a new 
multivariate cubature formulae Cjv (/) in the ball depending on a parameter 
G N which approximates the integral 

I f{x)dfiix) (1) 

Br 

for continuous functions / : Bf/ — )• C defined on the ball 

BR^lxeW^ : \x\ < R} , (2) 

where denotes the euclidean norm oi x = {xi, Xd) G . 

The exact definition of the polyharmonic interpolation formula and of the 
polyharmonic cubature formula Cn if) will be explained in Section [5] A major 
purpose of the present paper is to provide an error analysis for a class of functions 
on the ball Bji which exhibit a certain type of analytical behavior. 

Let us introduce the necessary notions and notations. Let 

S''-! := {x G M'' : \x\ = l} 

be the unit sphere endowed with the rotation invariant measure d9. We shall 
write a; e R"^ in spherical coordinates x ^ r6 with 6 € S'^~'^. Let Hk (K**) be 
the set of all harmonic homogeneous complex- valued polynomials of degree k. 
Then / S Hk (K'') is called a solid harmonic and the restriction of / to S"*"^ a 
spherical harmonic of degree k and we set 

flfe :=dim-Wfe(M''), (3) 

see [51], [in], [I], [H] for details. Throughout the paper we shall assume that 
the set of functions 

Ffe,, iM-^^C, for£-l,...,afe, (4) 
is an orthonormal basis of TLk {R'^) with respect to the scalar product 

Recall that due to the homogeneity of Y^.i [x) we have the identity Yk.i [x) — 
r^Yu (0) for x^r9. 

Our polyharmonic Interpolation and polyharmonic Cubature Cn (/) approx- 
imating the integral ([1} are based on the Laplace- Fourier series of the contin- 
uous function / : Bji — > C, defined by the formal expansion 

/(^^) = EE/MWi^MW (5) 

fe=0 1=1 
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where the Laplace-Fourier coefficient fkj (r) is defined by 

/fc/ (r)^ f f {re) Yk,i (6) dO (6) 

for any positive real number r with r < R and is defined in There 
is a strong interplay between algebraic and analytic properties of the function 
/ and those of the Laplace-Fourier coefficients fk,e- For example, if f (x) is a 
polynomial in the variable x = {xi, ...jX^) then the Laplace-Fourier coefficient 
fk,e is of the form fk/ (r) = r^pk,t (?'^) where pk^i is a univariate polynomial, 
see e.g. in [M] or [51]. Hence, the Laplace-Fourier series (O of a polynomial 
/ {x) is equal to 

dcg/ Qfc dcg/ ofc 

fc=0 £=1 fc=0 f=l 

where deg / is the total degree of / and pk,i is a univariate polynomial of degree 
< deg f — k. This representation is often called the Gauss representation. A 
similar formula is valid for a much larger class of functions. Let us recall that 
a function / ; G — > C defined on an open set G in M'* is called polyharmonic of 
order N ii f is 2N times continuously differentiable and 

A^u (x) = (8) 

for all X G G where A = -|- ... -I- is the Laplace operator and A the iV-th 
iterate of A. The theorem of Almansi states that for a polyharmonic function / 
of order N defined on the ball Bji = {x G M'' : < i?} there exist univariate 
polynomials pk,e (r) of degree < — 1 such that 

oc a/j oo ak 

f (^) = E E?'M(ki')>^,^ (^) = E E i^^i'pMd^i')^.,. w (9) 

fc=Of=l k=0 t=l 

where convergence of the sum is uniform on compact subsets of B^j, see e.g. 

EI], 0,0. 

To end this technical introduction, let us remind some estimates for spherical 
harmonics which we will need below: 

1. For every multiindex a and for every integer fc > 1 holds 

< C'fcl"l+'^ for6'e§'^-\ (10) 

see [in], p. 120. Here Dg is the multi-index notation for the derivative 
with respect to 9 ^ . 

2. A function / (0) defined on S''"^ is real analytic if its Laplace- Fourier ex- 
) satisfies 

rfcWIL,(s.-i) <C^e"'" forfc>0, 
for some constants G, > 0; see [5T| . 
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1.2 Complexification of the ball in MJ^, related to the ball 
of the Klein-Dirac quadric 

We want to study analytical extensions of functions / defined on the ball using 
the Laplace- Fourier series Our strategy is to require minimal assumptions 
on the functions /; thus instead of the standard approach where one works with 
functions / which are a priori analytically extendible to a fixed domain U in 
the complex space C'' (as in 2 ) we shall require only that we can extend the 
function x = rO i — > f [rO) to an analytic function z6 i — > f {z9) , so we only 
complexify the radial variable r to a complex variable z. Henceforth we will 
use the following terminological convention: If the function / {r9) possesses 
an analytic extension with respect to r we call the extended function / {z9) 
"r— analytic complexification" or "r— analytic continuation". Thus, for 
the r— complexification of the function / one should expect from equation ^ 
that the Laplace-Fourier coefficient fk.i (r) extends to an analytic function of one 
variable. Hence, we consider the analytically continued functions on domains 
in the set 

C X E'^-K 

We obtain the following important proposition. 
Proposition 1 The set of functions 

^'^'^'^-i k>0, e=l,2,...,ak, j = 0,l,...,iV-l / (11) 

is a basis for the multivariate polynomials which are polyharmonic of order N. 

The proof follows by representation ([7]). 

Remark 2 We will see later that the basis Bk.t-^N is a natural generalization of 
the basis {^"' }^q^ for the polynomials in the one- dimensional case. 

The main approach in the present paper is to consider the r— complexification 
of the functions / defined on the ball Bj^ C M.'^ in the form / {z9) . In partic- 
ular, for polyharmonic functions / we can provide a formal expression for the 
1 — complexification, using representation ([9|), by the following formula: 

k=o e=i k=o e=i 

The question of convergence will be addressed in the course of the paper. From 
this formula we make a crucial observation: the complexification / {z9) de- 
pends only on the values z9 and not on the coordinates of the pair (z, 0) . Indeed, 
this is due to the equality = {z9, z9) where for the vectors w,u ^ we have 

X — 

the non-Hermitian product by putting (u, w) := y ^ ujWj. Thus the function 
/ {z9) is defined on the space 

C X S'^-i/Zs = {z9 -.zeC, 9 e S'^"^} , 

where the factor in Z2 means identification of the points (z, 9) and (— z, —9) in 
C X S''"^. But the set C x §''~^/Z2 C is one of the possible representations 
of the famous Klein-Dirac quadric. 
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Definition 3 We define the Klein-Dirac quadric by putting 

KDQ := C X S'^-i/Za. (13) 

We shall see that the r— analytic continuation of the solutions of the poly- 
harmonic equation are in fact "the analytic functions" naturally defined on the 
Klein-Dirac quadric KDQ . The main interest of the present paper is devoted to 
the Function theory on the complexified ball Br in KDQ defined by 

Br := {z0 : \z\<R, 0e S''"^} = D^^ x S^^^VZa, (14) 

where D/j is the open disc of radius R in C, i.e. D^; {z G C : \z\ < R} ; for 
i? = 1 as usually one puts 

Di D and S = Bi. 

Some enlightening comments about the Klein-Dirac quadric KDQ defined 
in (|13p are in order. 

Remark 4 This quadric has been originally introduced in a special case by Felix 
Klein in his Erlangen program in 1870, where he put forth his correspondence 
between the lines in complex projective 3-space and a general quadric in pro- 
jective 5-space. The physical relevance of the quadric and the relation to the 
conformal motions of compactified Minkowski space-time had been exploited by 
Paul Dirac in 1936 J 13^ . The Klein-Dirac quadric plays an important role also 
in Twistor theory, where it is related to the complexified compactified Minkowski 
space 14^ - Apparently, the term "Klein-Dirac quadric" for arbitrary dimen- 
sion d, has been coined by the theoretical physicist I. Todorov, cf. e.g. 189^ . 
\4p]. In these references important aspects of the Function theory on the ball 
in KDQ have been considered in the context of Conformal Quantum Field The- 
ory (CFT). In the context of CFT Laurent expansions appear in a natural way 
as the field functions in the higher dimensional conformal vertex algebras ( us- 
ing a complex variable parametrization of compactified Minkowski space ); see in 
particular formula (4-. 43) in /^i/, as well as the references f39\/ . 140^ - 

Our main novelty will be a multivariate generalization of the classical Hardy 
space (D) called (polyharmonic) Hardy space on the ball Br to be in- 
troduced in Definition [TTl We will denote this space by {Br) . For simplicity 
sake we will restrict ourselves to considering the space (B) . 

Running ahead of the events, let us say that the name polyharmonic comes 
from the fact that {Br) may be obtained as a limit of the complexifications 
of the polyharmonic functions in the ball ([T^ : we take the closure of all finite 
sums of the type 

u{z,9) = Y,Y.''k,i(yZ^)z''Yk,i{e), (15) 
k=o e=i 

where Uk,i (•) are algebraic polynomials of degree < iV — 1; such functions u 
satisfy A^u {x) = 0. In view of the Polyharmonic Paradigm announced in |24j . 
the space {Br) generalizes the classical Hardy spaces which are obtained as 
limits of algebraic polynomials, where the degree of a polynomial is replaced by 
the degree of polyharmonicity. 
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The Hardy space (I^r) will be a Hilbert space and we will provide a 
Cauchy type kernel, which is the analog and a generalization to the Hua-Aronszajn 
kernel in the ball (cf. [5], p. 125, Corollary 1.1). Let us note that the last is a 
multidimensional generalization of the classical Cauchy kernel from Com- 
plex Analysis (more about Cauchy kernels see in [22], [32]). 

Remark 5 The reader familiar with the Cartan classification of classical do- 
mains, may remark that the boundary OBr is very close to the Shilov boundary 

{{e'^^re) : e [0,7r] , < r < i?, 6* G S'^-'^} 

of the so-called Cartan classical domain TZiv (called also "Lie-baW) equal 
to 

Bi -.^ {( + irj eC^ : ^,T] e R"^, q + iv) < l} where 




This has been considered from the point of view of several complex variables in 
the monograph of Hua lEJ^, and for the study of the polyharmonic functions of 
infinite order in the monographs (see in particular p. 59 and 126) and ]3^. 
The Hardy spaces defined in Definition \ll\ can be identified with the Hardy space 
of holomorphic functions on the Lie ball in C^, cf. 150]; this correspondence will 
be given a thorough consideration in J28^ . 

Remark 6 Let us define the annulus in the Klein-Dirac quadric KDQ as the 
set ^ 

Aa,b ■■= [zd e KDQ : a < |z| < 6, 0e S'^'^} . 

The Function theory on Aa^b would help to relate the present results to previous 
obtained by us. We have seen in \3Ui/ that an interesting, consistent and fruitful 
Function theory is available only on the set 

Aa,b^{{z,e):a<\z\<b, OeS''-'} 

which is a subset o/C x S'^^^. This is due to the fact that the r— analytic contin- 
uations of the solutions of the polyharmonic equations in the annulus Aa,b C M** 
live on the set Aa.b but not on the set Aa^b The point is that in the case of the 
annulus Aa.b we cannot identify the point z9 with {z,9) , in other words, {z,9) 
is not identified with (— z, —9) . 

The paper is organized as follows: in Section[2]we recall background material 
about the Hardy space H^ {V>r) . In Section [3] we introduce the polyharmonic 
Hardy space (Bb) on the ball B of the Klein-Dirac quadric. We prove that it 
is a Hilbert space, a maximum principle, and infinite-differentiability of the func- 
tions in H^ (Bb) . In Section [3T] we construct a Cauchy type kernel for {Bb) ■ 
In Section [5] we prove other main properties of the space ... which generalize 
similar properties of the one-dimensional Hardy spaces. In Section [5] we char- 
acterize the polyharmonic functions which are extendible to the Hardy space 

(Bb) ■ In Sectioniniwe prove some of the main results of the paper, about the 
error estimate of the polyharmonic interpolation, and about the polyharmonic 
interpolatory cubature formulas, generalizing the polyharmonic GauB-Jacobi 
cubature formulas introduced in |25) . [27) . 
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2 Classical Hardy spaces — a reminder 



Hardy spaces are a bridge between Harmonic and Complex Analysis. This is 
based on the fact that the Taylor coefBcients of a function / (x) on the real 
line K are at the same time the coefficients of an orthogonal expansion of the 
analytic continuation / (z) with respect to the basis {z^].^^ which is orthogonal 
on the circle. Thus in a certain sense the setting of the Hardy spaces represents 
a study of the properties of the real functions having Taylor coefficients with 
^ < oo by the methods of Complex Analysis. Since we are generalizing 
before all the Hardy space we will recall the main results about it. Let us 
put 

M{f;r) :^ i^l. J^^ \f {re^<f)\' 

Then the for every function / which is analytic in the disc D we define the 
Hardy space norm 

ll/llff2(D) := supAf (/;r) . 

Remark 7 For every analytic function f the function M (/; r) is an increasing 
function of r, cf. Theorem 17.6 in \44^ - 

A basic fact is that is a Hilbert space and may be identified with the 
limit values on the circle which coincide with (§^) . For every g & (§^) we 
have the norm defined by 




and the Fourier coefficients 

g (e*^) e-''''^dip for n e Z. 

Some of the main properties of the classical Hardy spaces are summarized 
in the following theorem (cf. Theorem 17.10 in [S]). 

Theorem 8 1. An analytic function f onJ} of the type 

oo 

= forzeB (16) 

3=0 

belongs to (D) if and only if 

oo 

in that case 

oo 
J=0 
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2. If f € (D) then f has radial limits f* (e"'') at almost all points on the 
circle S and f* G (§) . The Riesz condition holds, i.e. 

/;-0 for all J <0; (18) 
/; = /, for all > 0. 

The ^ approximation holds 

lim ^J^^\f i^^'") - f* (e'^) r = 0- 
The integral of Poisson and of Cauchy of f* recover f, i.e. 



2TTi Jr C ^ z 

where T is the positively oriented circle §. 

3. The mapping f i — y f* is an isometry of on the suhspace of L? (S) 
which consists of those g for which g{j) = for all j < 0. 

Let us recall the famous theorem of brothers F. and M. Riesz which concludes 
the absolute continuity of a Borel measure on §^ only from the annihilation of 
half of its Fourier coefficients, [HI Theorem 17.13]. 

Theorem 9 Let fi be a complex valued Borel measure on the circle §. // 

ivat 



d^i{t)^0 form>l, (19) 

then the measure fi is absolutely continuous with respect to the Lebesgue measure, 
i.e. the 
[0,2^]. 







I.e. there exists a function f* G (S) such that dix{t) — ^/* {t) dt for t G 



3 The polyharmonic Hardy space [B] on the 
ball of the Klein-Dirac quadric 

At first we observe that the basis functions hk.tj {r, 9) defined in (fTTjl have a 
natural r— analytic extension 

bk.Lj{z>^) ■■^z^^+'^Yk^eiO). (20) 

Let us present some heuristics for explaining our main goal: We want to 
define the Hardy space (ISr) as a space of functions which are uniform limits 
of sequences of complexified polynomials P {zO) on compacts of the ball Br in 
KDQ . For that reason we need an appropriate inner product. Let us make 
the important observation that there is a natural inner product where the 
basic functions (l20l) are orthogonal. Hence, for functions defined on the set 
§^ X W^^^ we introduce the inner product 

(/, 9).--=^ f r f (^'^' ^) JW^dipdO. (21) 
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The crux of our approach is the orthogonality of the basis functions bk^ij in 
PT|) on the boundary of the Klein-Dirac quadric §^ x §''~^/Z2, i.e. 

{bk,£,j,bk' ,e,j) — Sk,k'Se,e'Sjj' , (22) 

where the Kronecker symbol S means S^.p — 1 for a — (3 and for a ^ /3. 

This property is a remarkable generalization of the orthogonality of the basis 
on the circle §^ and traces the analogy to the one-dimensional case. 
Further, we provide some arguments about the proper definition of the norm 

of the prospective Hardy space. The objects of our polyharmonic Hardy space 
(B) will be functions / {z, 9) which are representable as infinite sums in the 
sense, and are absolutely and uniformly convergent on compacts with \z\ < 1: 



fc=0 £=1 \j=0 J 

= J2J2fkA^)YkA0) (24) 

k=0 1=1 

If fk,i,j = for all j > N and k < K ioi some ivT > 0, then the function / 
is a polynomial satisfying A^/ [x) — for all a; £ M''. Hence, formula ([25)) 
represents the r— complexification for the polynomials. 

In view of representation (1231) we introduce the following naturally defined 
subspaces of the classical Hardy spaces on the unit disc D C C: 

Definition 10 The "component spaces" H^'^ (B) C (B) consist of functions 
f e (B) having the representation f [z) — fi (z^) z'^. We put 

H^^^ (B) {/ (z) : / (z) = f, {z^) z\ f, G (B)} . (25) 

Thus the space H^''' (B) consists of the analytic functions in (B) having 
Taylor series 

oo 

/(z) = 5:a,z'=+^^ 

respectively the norm on 

(o) is the inherited from iJ^ (D) . 

Definition 11 We define the polyharmonic Hardy space (B) on the unit 
ball B = Bi of the Klein-Dirac quadric defined in |_?.^[ ), as the space of func- 
tions f given by the Laplace- Fourier series fig) with coefficients fk,e G H'^''' (B) 
satisfying 



11/11 



\ EE ii/fc.^ii^^(D) < °°- 

\ k=o e=i 



Remark 12 The reader may note that Definition \11\ mimics the definition of 
the classical Hardy spaces which are obtained as the closure of the polynomials. 

Now let / e (B) . By Definition [III since aU fk,£ G H'^ (B) it follows that 
for r -J. 1- and z = re*"^ aU fk,i {re"<') have limiting values f^/, {e'"^) in L2 (S^) , 
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hence 

oo ak 
k=Q i=l 

This implies 



k=0 1=1 



p2n 








/ \n 




Jo 





(26) 



^lim| j \f{re^'^,9)\^d^de = Y,Y^J\fl,{e^'^)\'d^<oo, 

Sd-i '^=0 

hence, it follows that / G L2 (B x S'^"!) , and due to also f £ L2 (B) . 

Note that as in the classical Hardy spaces the inner product ((2T|) is good 
only for the polynomials but might not be well defined for arbitrary functions 
having bad boundary behavior. For that reason we have to change the definition 
of this inner product. 



Definition 13 We put 



if, 9) 



1 



S<i-1 Jo 



2-77 



f {re"^,e) g{re^v,e)dipd0. (27) 



The following theorem justifies our arguments above and is an analog to 
results for the classical Hardy spaces, cf. Theorem[5](or Theorem 17.10 in 1441 ) . 

Theorem 14 1. The space [B) is complete. 

2. It coincides with the space of functions f having representation 



00 afc 00 



/(z,0) = EEE/M;.^'+'-''i^M (0) 

k=a i=i j=o 



with coefficients satisfying 



00 afc 00 



EEEi/M.i^<- 

fc=0 1=1 j=0 



(28) 



(29) 



3. The norm of f is given by 



00 ttfc 00 



(eeeia^^.1 

\^k=0 1=1 j=0 



1/2 



The following equality holds, for z — re^"^ , 



k=Ol=l 



2j 



3=0 



00 ak, OQ 



^■^^EEEiA-.^j 

fc=o 1=1 j=n 



(30) 

4-. Every element f £ (B) is the limit of a sequence of polynomials 
Pn G V which satisfy 

A^Pn (x) = forxeBc ISJ^. 
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We will prove only some of the above statements. 
Proof. 1. It is obvious by a direct estimation that condition / S [B) implies 
the absolute convergence of the series ([25]) on every compact if C B x S'^""'^. 
Hence, every function defined by the series (051) and satisfying also satisfies 



c» Qfc oo oo oo 



E E ii^,^ii^^(B.) = H E E E i/M.f r^''^'^' =j:ej: i/^,^;/ < 

k=o e=i k=o i=i j=a k=o e=i j=o 

which shows that the space (B) is a tensor sum of the spaces {/fc,«:j}j>Q 
taken over all indices (fc, i) . Now, in a standard manner, the completeness of 
the space (B) follows from the completeness of the space of sequences £'^. 

4. We have to see that every function defined by the series (pS)) and satisfying 
([2^ is a limit in the norm of (B) of a sequence of polynomials. Indeed, for 
every £ > we may find integers fci > and > such that 

oo dk oc 

E E E i/m;/<- 

fc=fci+l i=l j=N+l 

Then the polynomial P defined by 

fci dk N 

p{x) :=EEE/m;.'''+'^>^,^ ie) 

k=0 1=1 j=o 

satisfies with z — re*''' the following: 

27r „ ^. „ 



- / \f{ze)-p{ze)U^de = E E E l/- 

^ 1 1, I 1 fi 1 A AT I 1 



gd-i k=ki+ie=ij=N+i 

OO C^fc oo 

< E E E i/M.i^<- 

k=ki + l 1=1 i=N+l 

This implies 

||/-P||«.(i5)<£, 

which ends the proof. ■ 

Remark 15 The. essence of Theorem \14\ is that, as in the classical Hardy 
spaces, only the information about f in the real domain, provided by the Laplace- 
Fourier coefficients 

/ (^) = E E A-.^ -"^^^^ (^) - E E E hA^r'^^'^Y'^-^ (0) ^ 

k=0 1=1 k=0 1=1 j=Q 

determine when does f belong to (B) . All we need to know is that they satisfy 
the convergence condition i29\) . J2'k'=oJ2'e=i^'jLo \ fk,i,j\ < oo. 

Theorem [Ml suggests the representation 

k,£ 
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where the tensor sum is understood in the sense of equahty ([M)) . and the prime 
in the symbol ^ means that only sums are taken which are convergent as 

The following result follows from the representation in (PS)) . 



Proposition 16 For every function f G H^ '^ (D) the following formula of 
Cauchy type holds: 



here C e ID). 



^^f(z)dz; 



3.1 Cauchy type kernel for H'^ (B) and Hua-Aronszajn type 
formula 

The orthogonality of the basis {bkj.j} hints us to construct a Cauchy type 
kernel and a corresponding formula which reproduces the multivariate poly- 
nomials by using their values on the set §^ x §'^"^/Z2. By exploiting the above 
orthogonality, we obtain such formula easily, following the general principles of 
constructing kernels, by putting: 

OO Qfc 

j>0 k=0 t=l 

OO Ofc 

= E E E c'^'^'^fe,^ (^') ^''^"Y^.^ (^) ■ 

j>0 k=Q 1=1 

Note that this kernel is absolutely convergent for every ^ and 9 with |C| < 1, 
9' E and z — e*''', due to the estimates for the spherical harmonics (fTU]) . By 

the orthogonality property (j22p . for every polynomial P we obtain the Cauchy 
type formula formula 

P{Ce')-{K{CO';-),P{-))^ (32) 
1 



27r 



S'i-i Jo 



K {C,9';z,9) P{z,9)dipd9. 



Let us remark that no such formula is available in the real domain. Formula 
(15^ is a strong motivation to consider further the consequences of the inner 
product (1^ . 
Obviously, 

OO afc 

K (C, 9'; z,9)^J2 (C^)'' E E (C^)' ^fe.^ (^') ^^.^ W 

j>0 fc=0 £=1 

, OO Ofc 

= TT7^EE(C-)'^m(^')^m(^) 

which shows the absolute convergence for |Cz| < 1. Let us recall that the usual 
Poisson kernel (see ^54j, chapter 2, Theorem 1.9) is given by 

1—2 OO dk 

^' I. ~[nd = E E ^'^fc^^ (^) ^fc.^ (^') < 1 

b'-T-fc' fe=0£=l 
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in every dimension d > 2. This expression is obviously close to the above expres- 
sion for K (C, 9'; z, 6^ . We will apply the idea for the r— complexification to the 
kernel Kp (r, 9, 9') and we will relate it to the Cauchy type kernel K (C, 6'; z, 0) . 

Proposition 17 For every complex number w with \w\ < 1 the following equal- 
ity holds: 

oo oo 1 2 

y: e ^'^^^^ (^) ^^^^ (^') - E ^'4'^ (^') = — — (33) 

fe=Ofci fc=o (l-2u;(0,^?')+ii;2)^ 

where zf^ (9'^ are the zonal harmonics, see J54l- Let us put cos0 — {9,9'") . 
For d = 2 we have 

( n ni\ 1 1 X"^ u COS k(j) 1 l~w^ 

Kp(w,9,9')^ h-> w ^ = TT (34) 

^ 27r TT ^ TT 27r 1 - 2w cos + w2 ^ ^ 

/c— 1 

and for d > 2, 

ifp {w, 9, 9') = E ^'cfe,d^fc (cos 0) (35) 

fe=0 

where 

I 2k + d-2 d-2 

Ck,d 1 — ;^ — , A — — — 

LOd-i d — Z 2 

and are the Legendre polynomials. 

Hence, the r—complexification of the Poisson kernel Kp (r,9,9'^ is given 

by 

1 — 7/;^ 

Kp {w, 9, 9') = (36) 

{l-2w{9,9') +w^)' 



1 — w 



2 



{{1 - e"i'w) {I - e-"f'w)f^^' 

The Cauchy type kernel K (C, 9'; z, 9^ defined in 1131]) is related to the r—complexification 
of the Poisson kernel by the equality 

K{C,9'-z,9)^ \ Kp{Cz,9,9') (37) 

1 — 

1 



{l~2Cz{9,9')+ez^y 

Proof. The proof uses the following results from [54] . 

1. Formula ([M)) follows from formula (2.7) in chapter 4 in |54| . 

2. Formula psp follows from Lemma 2.8, Theorem 2.10 and Theorem 2.14 
in chapter 4 in [54 . 

By the estimates for the spherical harmonics in (jlOD it follows that the 
left-hand side of (j33p is absolutely and uniformly convergent on every compact 
K cBx S"*"^ X S'^"^ Representation ^ follows from formula 



1 - 2w cos = (1 - e"l'w) (1 - e-"l'w) . (38) 
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The definition of tfie kernel K (C, 9' \ z, 6*) shows that for all polynomials P 
with real coefficients holds 



1 



S<i-1 JO 



2tt 



P(C6I') = ^;^/ / K {C,e';z,9) Piz,e)dipde forz = e^^. 



Since z — z ^, a change in the integration ip ^ —ip shows that 
PW^T^I I K{C,e';z,e)P~{^dp,d9 
' ' K [C,9';-,9] P{z,9)d<pd9 



S'i-i Jo 



— f I -K ((:,9';-,9] P{z,9)dzd9 
lujd J%d-i Jr^ z \ z J 



2Triuj 

where Fi is the positively oriented circle S^. On the other hand we see that 



-KIC, 



{z^-2Qz{9,9')^CY 

The last by definition is up to a factor the Hua-Aronszajn kernel, see [2], p. 126. 

The above identities motivate the definition of the well-known Hua-Aronszajn 
kernel H (^(,9';z,9) given by 

H{C,9';z,9)^— — (39) 

{e~2Cz{9,9')+z^y 

where uJd = n'^^^/T {d/2) is the surface of the sphere (see [2], p. 122 Theorem 
1.1, and p. 126, Remark 1.4, where up to a factor it is called Cauchy kernel for 
the Cartan classical domain TZjy)- 

From above we see that the following equality holds: 

H{C,9';z,9)^--k(c,9';-,9)=- ] £ (C/zf Y,,, {9') Y^j (9) 

(40) 

1 1 1 



'^'^^(l-2C/z(0,6'') + cV^^)' 



4 Main properties of the Hardy spaces [B) 

In the next theorem we provide a generalization of the classical boundary value 
properties of the Hardy spaces if^, see e.g. Theorem[5](or Theorem 17.10, 17.12 
and 17.13 in [H]). 
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Theorem 18 Let f e H'^ (B) . 

1. Fatou type theorem: For r — > 1^, and for almost all (p G [0,27r] 
and 9 e S''^^, the function f {re^'^O) has a radial limit which we denote by 
f* {e"^, 9) , and which satisfies f* {e"^, 9) e L'^ {S x S'^-^) . If the expansion of 
the function f is given by 



oo a/^ oo 



/(z,0)=^^^/M;,z'=+'-'>fc,£W for all \z\<l,9eS''-\ 
k=o 1=1 j=a 

then f* (e*'^,^^) is given by the Laplace- Fourier series 

oo afc oo 

/* (e^^ 0) = E E E (^) ^ = e^^. (41) 



fc=0 i=l j=0 

2. The following limiting relation holds, 



lim— / / \f (re''^9) - f* (e''<'9)\^ d(pd9 ^0. 

r^l 27r Jq Jgd-l 

3. Let the Laplace- Fourier series of the function f* (e*'^, 9) G (§ x §^~^) 
be given by 

oo ttfc oo afc oo 

/* (e^^ ^) = E E /fc/ (e^^) ^^.^ (^) = E E E (^) V ^ = 

/c=0 £=1 fc=0 ^=1 j=0 

Then the Fourier coefficients of the functions f^ £ {^^^) satisfy the following zero 
conditions, which we call Riesz type conditions: 

f*kA3=^ forallj^k,k + 2,k + 4,.... (42) 

These conditions are equivalent to (the usual form of Riesz conditions) 

271- p 

/ f*{e'f,9)e-'^'^Yk,e{9)dipd9 = for all j k, k + 2, k + 4, ... . 

JS'i-i 

(43) 

4. The following Cauchy-Hua-Aronszajn formula holds 



f{C,9') = ^ [ [ H (C, 9'- z, 9) f* (z, 9) dzde 



(44) 



^ ^ k(c,9';-,9\ f*{z,9)d^d9 



"'S' 



^ ^ Pr2{2ip~2Lp') [ Kp (-,9,9' ] f* {z,9)d9dip, 



2TTU}d Jo 



where we use the notations z = e^^ and = re^"^ , the kernel K is the Cauchy 
type kernel LSI]) , and Pr {(p) is the usual two-dimensional Poisson kernel. We 
call the kernel 

Pr2{2(p-2ip')Kp(-,9,9'] (45) 
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the modified Poisson type kernel. Here the contour Ti is the positively 
oriented circle or a bigger contour which encircles it. 

5. Dirichlet problem with L'^ data: If a function f* (e*'^,6') G (S x S'^^^) 
satisfies the Riesz type conditions J^^l j then there exists an unique function 
f e H-^ (B) which has as a "non-tangential limit" the function f* (e"^,0) in the 
sense ^ 

lim / / \f (re''^,9) ~ f* (e''^,e)f dipde = 0. 
— >1 JS'i-i Jo 

6. Every function f* £ x S''^^) which satisfies the Riesz tyep con- 
ditions l^iip belongs to the space (S^ x S"^'^^ /'Z2) ■ The map f — 5- /* is an 
isometry between (B) and the subspace of (E>^ x §''^^/Z2) . 

7. The norm is given by 



1 



27r 



Proof. 1. We will omit the proof of this non-trivial but standardly proved 
statement. 

2. By Theorem [Til the function / is represented by a series (|28p satisfying 
Sfc^oSfci X^jlo IfkAjl^ < Let us consider the function 

00 ak 00 

/* (e^^ 0) ■■=J2J2T. hAjz'^^^'Yk,, (9) for z^e''^. 

k=0 1=1 j=0 

Since the functions {e*-''^Yfc ^ (6) : j G Z, all {k,£)} form a basis of the space 
L2 (§1 x S'^-i) , it follows that /* G (§1 x S"^-^) . By Parseval's theorem we 
see that 

•^^ k=o e=i j=o 

It is easy to see that for p — > 1~ the last tends to 0. 

The proof of 3.) is evident due to the expansion (j4ip . 

4. From formula gO]) it follows that for every function g (e'*^, O) G x S''"^) 

holds 

' f f ^'-'f*^^^%zde = 



K {C,0';z,0) f* {z,0)dipd9 for 2 = 6'"^. 



27ru;„ 



d Jo JS''-i 



From the series representation of / and /*, by the definition of the kernel 
K{C,e';z,0) it follows 

^T'' Jo JS<*-1 

Further, since 

00 

P,. - = ^ rljle*j('P'-¥') 



j = -oo 



16 



and 



the product {'2(p — 2ip') Kp (^^,9,9'^ has terms which in the integral will 
vanish due to the Riesz conditions satisfied by /*, and there will remain only 
the terms corresponding to K (C, 9'; z, 9) , which proves the formula. 

5. Every function /* (e*'^,0) which satisfies the Riesz type conditions (jH]) 
has a Laplace- Fourier series of the type (|4ip . say 

oo afc oo 

/* (e^^ ^) = E E E ^kAj^'^''Yk,i {9) for z = e''^. 
fc=o e=i j=o 

Since /* {e''^,9) e (S^ x S'^-^) , the coefiicients satisfy J^T^o T,Zi T,JLo |afe,^;i 
oo, hence, we may define the function 

oo ak oo 
k=0 1=1 j=0 

the limiting property follows as in item 2). 

Finally, item 6.) follows from the above arguments. 
■ 

Theorem llSI shows that we may solve Boundary Value Problems in the spaces 
(B) which is an essential advantage over the situation with the holomorphic 

functions in and alternative definitions of Hardy space in several dimensions, 

see [Sg, 05], [51], gS]. 

By Theorem [T8l the following space: 

{/* G (§ X S'^^i) : /* satisfies the Riesz type condition (gH)} 

is isomorphic to the space (B) . The Cauchy-Hua-Aronszajn formula pi)) 
generalizes the Cauchy formula in C and the Poisson formula in R'' at the same 
time. 

Remark 19 Let us formulate a conjecture about an analog to brothers' Riesz 
theorem: Let the complex valued Borel measure ji {ip, 9) be given on E>^ x S'^~^ 
with ip e [0,27r] and 9 G S"^^^. Assume that for all indices {k,£) holds 

z^Yk,i{9)dfi{ip,9) =0 for j ^ k,k + 2,k + 4,... (46) 

Ls the measure fi "absolutely continuous", i.e. does there exist a function f* 
which is in Li (S^ x S'^^-'^) such that d/i {(p, 9) = f* {ip, 9) d(pd9 ? It seems that 
the answer in this form is negative, but a positive answer needs some addi- 
tional properties of the measure pi. A thorough discussion to this question will 
be considered in \2Sj . Let us remark that a genuine analog of the brothers Riesz 
theorem is difficult to achieve for all approaches to Hardy spaces, cf. J19f . 

EI, m, 
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Remark 20 Another conjecture about the boundary properties of the remark- 
able Poisson type kernel in |^^[ ) generalizing a classical situation is the following: 
Assume that the function g (e^,0) belongs to C (S^ x S'^~^/'Ii2) o-nd satisfies the 
brothers Riesz conditions {4^ - Then by means of formula (or equivalently, 
by {4-4^ ) we may define a function Fg {z,9) on the interior (for |z| = r < 1 ) 
of the ball B of the Klein-Dirac quadric. Let us put Fg (e*''',^) = g {e^f ,9^ for 
r = 1. We conjecture that the function Fg is continuos on the closure B. Let 
us note that in the classical case the Poisson kernel is used to prove similar 
statement, see chapter 2, Theorem 1.9 in ,'5^/. Here we expect that the modified 
Poisson kernel \45^ will of central importance for the solution. 



4.1 Maximum principle 

In the classical case of the Hardy spaces, the niaximum principle is intimately 
related to the Cauchy formula in C or to the Poisson formula in R'^ (see the 
proof of the completeness of H'p in Remark 17.8 in [H]). A weak form of 
maximum principle alows to prove that the elements of are uniform limits 
of polynomials on compact subsets of D. Here we prove analog to this for the 
polyharmonic Hardy space (B) . In the next theorem we see that the explicit 
form for the Cauchy-Hua-Aronszajn kernel is essential for proving a maximum 
principle. 

Theorem 21 Let f e (B) . For every q with < q < 1 we have the following 
(weak) maximum principle holds: 

1/ m\ < (1 - q)-' UWh^b) for all Id <q, 9e 

More generally, for every mixed derivative D" g with respect to the variables C 
and 9, we have the maximum principle 



\D"f{CO)\<C, X (l-g)-'^-l"l 
for all Id < q, and 9 G S'^-^; 



here the constant Ci > is independent of a. Respectively, for real C — ^ with 
X = r9 this gives an estimate for D"f (x) . 

Proof. Indeed, for |d < q and z = e'*^, if we put w = (/z, by the Hua-Aronszajn 
formula (l44l). it follows 



\fm\<^ I 



27r 



1 



K[(:,9';- 



\f* {z9') \ \dz\d9' 



< 



1 

2^ 
1 

2^ 



S'i-i Jo 



(l - 2^(61,61') 

dipd9' 



1/2 



dipd9' 



1/2 



Here we apply (j38p with cost/j = (^9, 9'^ , which implies the inequality 



2w{9,9')+w^) \ = 1(1 



s^'^w) (1 



e~^^w)\>{\-q)\ 
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hence, 



(1 - 2w {6, e') +w^)^ =1(1- e''^^) (1 - e-^'I'w) I 



>(i-'zr 



By Theorem [TH] this imphes 



1/ (Ce)l <(!-'?) " 



which ends the proof of the first part of our statement. 

In the same way, by differentiating under the sign of the integral in the 
Hua-Aronszajn formula, we obtain the estimate for the derivatives, 




X 



By the Leibniz differentiation formula, the differentiation with respect to C and 
6 gives the following estimate 



< 



< C 



d+\c 



- - 1 



X \a\ 



and this ends the proof. 
■ 

We have the following immediate corollary about the regularity of the func- 
tions in the space {B) . 

Corollary 22 The functions in H'^ (B) belong to C°° (B) . 



The proof follows from the maximum principle in Theorem[2T] since every / e 
(B) and the derivatives of / are uniform limits of a sequence of polynomials 

Pn {z9) and the respecitve derivatives of Pat {z9) on every compact sets A'xS''^^ 

where the compact K dJi. 



4.2 Real analytic functions and the space (B) 

As is the case with the classical analytic functions the elements of (B) are 
real analytic which we prove in the next theorem. 

Theorem 23 Let F G (B) . Then the function f [x) = F {r9) is real- analytic 
in the hall B dW^. 

Proof. We use the maximum principle: By a change of the variables, and a 
routine estimation with the Leibnitz rule, we obtain 

\D'^f(x)\ < C max 
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The last is estimated by the maximum principle on compacts, by Theorem 1211 
which gives 



1^1 



for r < q < 1. 



The statement of the theorem follows by an equivalent definition of real-analytic 
function (cf. 0, p. 17). 



5 BVPs for the polyharmonic operator and 
the spaces (B) 

5.1 The one-dimensional case 

Let us provide some hcuristical observations from the one-dimensional case. 

As we already said, the Hardy space setting provides interpretation to real 
Taylor series / (t) = J2jLo o-j^^ for which we have Yl'jLo 1^ < oo as the Fourier 
series X^jlo o.j^^''^ in the Hardy space i?^ (D) . Thus the Hardy spaces provide 
a different way to encode the information of real-domain data. There is an 
alternative way to encode the information, by means of boundary data, e.g. by 
means of 

It is more convenient to consider further differential operators defined on the 
functions on §. 

Let be a function defined on S. Indeed, if we have a function F which is 
representable as a series 



oo 

ij'-P 



3=0 



then we obtain 



where we have put 



i{j-l)ip 



We want to investigate for which sequences {cj}°^Q and {rfj}°°^Q there exists 
a function F e H'^ (D) such that 

D^^F{-l)^cj, D'^^F{l) = dj for j = 0,1,2,... . 

since D^^ is the operator A-' in the one-dimensional case. The matrix of the 
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system for the unknown coefficients aj is given by 
/ 1 1 



1 -1 



111 


1 1 


1 -1 1 


-1 1 


2-1 3-2 4-3 


5-4 6-5 


2-1 -3-2 4-3 


-5-4 6-5 


4! 


5! 6!/2! 


4! 


-5! 6!/2! 



V 

We see that even in the one-dimensional case the map /* 
trivial; this is in fact the map to the Taylor coefficients 



A-^/igd-i is not 



oo 

{D'^F (-1) , D'^F (1)}^.^^ ^ {f(^^ (0)} ^F{z) = Y^ 



J^(J) (0) J 



j=0 

5.2 Which polyharmonic functions are in (B) ? 

The polyharmonic function in the one-dimension case satisfy d^P/v-i (i) /dt^ — 
and are polynomials. Their complexifications P/v {z) belong to the Hardy 
spaces in the disc. On the other hand, not all polyharmonic function (of fixed 
finite order N) belong to the polyharmonic Hardy spaces (B) . In the present 
section we characterize those polyharmonic functions which belong to (S) . 

We will characterize the polyharmonic functions by means of their boundary 
properties. The main point is that the polyharmonic functions in a domain 
D C Mf^ (and more general domains) may be " parametrized" by considering the 
Dirichlet problem 





9j 



in D 
on dD 



for j = 0,1,2,..., iV-1. 



By means of the classical Green formulas we may find u (x) for x E D. On the 
other hand, if the function u which has a r— complexification to ball B of the 
Klein-Dirac quadric, then by means of Cauchy type formula (j44p we may recover 
u {z, 9) using its values u {e^'^ , 0) on the boundary dB = S x §''~^/Z2. Thus we 
see the boundary values u{e^'^,9^ provide an alternative parametrization for 
the polyharmonic functions. It is essential, as in the one-dimensional case, to 
provide a relation between the boundary data {ffjl^o^ ^^'^ (^''''' ^) • 

The following theorem establishes the link between the polyharmonic EVPs 
and the boundary values of the elements in (B) . 

Theorem 24 Let u be a polyharmonic function of order N > 1 in the ball 
B C M'', i.e. A^u{x) — in B. Then the r— complexification of u satisfies 
u G (B) if and only if 



1, 



for j = 0,1,..., N 

where (S'^^^) denotes the Sobolev space of exponent s. 

Proof. First recall some properties of the polyharmonic operator A^. 
following operator is defined by 



The 



L(k)f 



1 d 



r 



d+2fc-l 



dr 



If 

k 
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see p. 152 in [53]. For / (r) = r'^+^'', k,s >0, we obtain 

(47) 

Hence, for every j < s holds 
Hk)[^'^''] = 4^s(s-l)---(s-j + l)x (48) 

xg + fc + .-l)...(^^ + fc + .-j) 

Let us put 
Obviously, 

7^;, = for s<j- 1. 
Let the function u (x) have the expansion 

oo ttfc A^— 1 
fe=0 £=1 j=0 

and for the boundary data we have the Laplace-Fourier series expansion 

CXD Ofc 

fc=0 i=l 

Hence, by formula in [24], p. 165 for the Laplace operator, we have 

oo ak N — l 

(6) = A^u (a:) 1,^1 = E E E ^f^^^-^^LTk) [r"^''] YkJ (0) 

k^O £^1 j^O 

oo N—1 

= EEE"M;.7l,>^M 

fe=0 i=l j=0 

By comparing the coefficients we obtain the system of equations 

5m = E for m = 0, 1, iV - 1. (50) 

j=m 

Hence, the map u -f-s- {<7j}^q^ is given by the infinitely many matrices for 
A; = 0,1, 2,..., 





( 


1 


1 


1 


1 


\ 









7^2 • 














72,2 ■ 








[ 








■ 


■ 7Af-l,Af- 


1 / 



^d+2k- 



dr 



1 d 



k-i± 

dr 



^d+2fc-lJl r 2sl 
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The result follows from the Lemma below, and by application of the 
Cauchy-Bunyakowski-Schwarz inequality to the norm 



N — l oo 
j=0 fc=0 1=1 



with equations 



We have the following technical lemma. 
Lemma 25 The solution of the system i5U\) is given by 

X 



=] 



where Dk — detC/^, and -D^!'"' are the appropriate minors of Uk given by the 
Cramer's rule. These coefficients satisfy 

> Cpj for fc oo, 

for appropriate constants Cpj. 

Proof. In the system ((50)) we divide the j— th row by (f + kY . On the 
right-hand side we obtain the vector , ^ ,■ , and the coefficients of the 

system of the left-hand side (|50)) become equal to 

. for 1 < s < TV - 1. 

ii + kY 

From formula (|48|) for the constants 'jj ^ we see that 



Fj_s for k oo, 



ii + kY 

for some appropriate constants Fj^g which are non-zero. This ends the proof 
since the minors D,'^ and Di. are computed through the values , ^ . 
■ 

We see that TheoremlMlprovides us with a large class of functions defined on 
the ball B cR'^ which are extendable to the ball B of the Klein-Dirac quadric. 
It is also possible to consider functions which are in a certain sense polyharmonic 
of infinite order, as those studied by Aronszajn, Lelong, Avanissian, and others, 
cf. the references in [2], [3], [29] . 

6 Error estimate of Polyharmonic Interpolation 
and Cubature formulas 

The topic of estimation of quadrature formulas for analytic functions is a widely 
studied one. Beyond the classical monographs |35j . |12| . we provide further and 
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more recent publications, as [S], [IS], [T7], [13], [SI]) [31] ■ No references may 
be found though for the multivariate case, for cubature formulas, even in the 
fundamental monographs as [51], [5S], [55]; see also the recent survey [TUj . 

Our main framework of Interpolation and Cubature was defined in [57] , [IH] , 
[25j . It has further brought to life the multivariate complexification and the 
polyharmonic Hardy spaces. 

We will consider first polyharmonic interpolation, and as second, polyhar- 
monic Cubature formulas for approximating integrals of the type 

/ f{x)dii{x) 
J B 

over the unit ball B C M'', where ^ {x) is a signed measure of special type. 
The main purpose of the present section is to find error estimates for the In- 
terpolation and Cubature in the case of functions / G [B) . Interpolation 
and Quadrature are intimately related in the one-dimensional case, and we will 
demonstrate similar relation in the multivariate case in our setting. 



6.1 One-dimensional case 

First, following the classical scheme outlined in [T^], p. 303 — 306, (see also 
chapter 12 in [32]), we will remind how one finds the error for the classical 
Quadrature formulas. 

Let the points to, ti, belong to the interval [a, h] . We define 

{z) = {z - to) (z - ti) ■ ■ ■ {z - tN) . 

Let the function / be analytic in a simply connected (open) domain D d C 
containing the interval [a, b] with boundary dD = T. Then the interpolation 
polynomial P/v {t) satisfying P/v {tj) ~ f {tj) for j = 0, 1, is given by 

j-' uj' (tj) (z ~ tj) 2Tn J oj{t){t-z) 

where V is considered as a contour oriented counterclockwise. Hence, the re- 
mainder is 

2Tn J u) (t) (t — z) 
r 

Now, if /X is a non-negative Stieltjes measure, say d/i (t) = w (t) dt, the quadra- 
ture formula 



J /(i)d/.(i)«^A,/(t,) 



is called interpolatory quadrature formulas of degree N if it satisfies the 
following equality 

N 



Y,X,Q{t,)^ I PNit)df,{t) 
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for every polynomial Qn oi degree < A^. This implies that 

of. [12, P- 303, or Krylov, [3S], chapter 12. Hence, for the error of such formula 
we obtain 

E if) / (/ (.) - P. (.)) d^i^) = hj{j ^wi^A • 

a a \V / 

This may be directly estimated by 

b 
a 

where d := min-,- (dist (tj, T)) , D := maxj (dist (tj^ {a, b})) , S-p := mindist ([a, h] , V) 
and Lr is the length of the contour F. 

6.2 Multivariate Interpolation 

Now we consider the multivariate case. 

Let us assume that the number b satisfies 

< 6 < 1. 

Choose the domain D = Bi and a function / G (B) , assuming that / has 
the expansion 

fe=0 1=1 

Let > be a fixed integer. We will consider polyharmonic interpolation 
which has been studied in |18j . Let the points {fk,e;j}^^Q belong to the interval 
[0, b] . We consider the following series: 

PN(z,e)^J2J2 Pfc.^ (^') ^'^^",1 {9) , (51) 

i; = 1 = 1 

where and pk,i are polynomials of degree < N satisfying the interpolation con- 
ditions 

Pk,t {rlfy) = fk,i {rli;j) for j = 0, 1, N. 

We prove below that the series ([ST|) is convergent. 
For the remainder of this interpolation we have 

oo afc 

/ (z, 0) - Pn {z, 0)=Y,T. [/'^^^ (^') - P^'^ ^'^^.^ (^) ■ 
fc=o e=i 

Now define as above the functions 

UJkJ (z) = {z- rl i.a) {z - rl^.^) ■ ■ ■ {z - r^^.^) , 
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and consider the oriented contour F (t) — e'* for t 6 [o, 2tt] . 
For all z with \z\ < b and 9 G S'^^^, and obtain the estimate 



\fiz,e)-PNiz,i 



< 



< 



< 



< 



(z^) //c,f (r^) 2t 



CO Ofc 



< 



(52) 



EE 

k=0 l=\ 
CO 

EE 

fe=0 1=1 

CO Ofc 

EE 



< C 



< c 



r 

2^+2 



^kA{T^)iT^-Z^) 

iOk,i (^') Ikj (t^) 2 



dr 



(1-6) 



Ar+2 



(1 - b)^+^ 
V-^ A;=0 e=l 



EE 1/^, 



Idrl X 6'=|yM {6)1 
(r2)| |dr| X 



2^+2 



(1 - 6)^+' 



■EE 11/^. 



Since / G i?^ (B) the last inequality shows, after application of Cauchy-Bunyakovski- 
Schwarz inequality, that the series J2kLo J2e=i [fk,e (z^) - Pk,e {z'^)] z'^Yk,t (0) 
is absolutely and uniformly convergent. Hence, the series (|51l) representing the 
function P]\f {z,9) is also such. 

Above we have outlined the most important arguments for proving the fol- 
lowing: 



Theorem 26 Let f G (B) . Let the points {fk,e;j}^^Q belong to the interval 
[0,&], where b < 1. Then the function Pf^{z,9) defined by the series i51\) is 
polyharmonic of order in the ball C ffi.'' and belongs to the polyharmonic 

Hardy space (Bb) , while the following inequality holds: 



l|-P'v||//2(B^) < CN,b ll/lli/2 



(B) 



6.3 Multivariate Polyharmonic Cubature 

The class of pseudo-positive measures used for our cubature formula Cat (/) 
is now defined in the following way: a signed measure /i with support in B^i C 
M'^ is pseudo-positive with respect to the orthonormal basis Yk,£,£ — l,...,a/c, 
fc G No if the inequality 

/ h{\x\)Yk,eix)d^iix)>0 (53) 

holds for every non-negative continuous function h : [a, b] — >■ [0, oo) and for 
all k G No, i = 1,2, ...,afe. Let us note that every signed measure d^ with 
bounded variation may be represented (non-uniquely) as a difference of two 
pseudo-positive measures. We refer to |27j for instructive examples of pseudo- 
positive measures. 
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Let the pseudo-positive (signed) measure be given in the ball Bi, C R'^. 
For all indices {k,£) the component measures are defined by 

df^k e (r) ■■= I Yk.i {9) dfi {re) > for aU r e [0, b] ; (54) 



here the integral is symbolical with respect to the variables 9. Rigorously, the 
component measure c?/ij, ^ (?') is defined for the functions g (r) on the interval 
[0, 6] by means of the equality 

b 

J 9 (r) diik,e (r) J 9 (r) Yk,i (9) d^x (x) ; 

cf. m. m- 

In [3S], [37], we have considered a special type of Cubature formula, the 
so-called polyharmonic Gauss-Jacobi Cubature formula. Here however we will 
consider more generally, interpolatory polyharmonic Cubature formulas and will 
prove their convergence and error estimate for them. The case of the annulus 
has been considered by us in [30] . 

Let us fix {k,£). We assume that there exist points tk^e-j, j — 0,1, ...,A^, 
belonging to the interval [0, b] , and numbers {^k,e;j}j^i , such that the following 
interpolatory quadrature formula holds: 

^. N 

/ Q (r) d^j,^^ (r) = ^ Xkj;jQ {tk,i;j) for every Q £ Vk,N; (55) 
here the set Vk^N is given by 

We define the polyharmonic interpolatory Cubature formula by 

OO Qfc N 

if) := E E E • (56) 

k=0 i=l j=0 

For the interpolation polyharmonic function Pn defined in (|5ip we obtain equal- 
ity 

Pn {x) dfi (x) = Cn {Pn) ■ 



Hence, the remainder of the polyharmonic Cubature formula is given by 
E{f) ^ [ {f{r,9)~PN{r,9))dfi{z,9) 



EE/ [hAr')-P'^Ar')]r'd^x,^iir), 



k=0 1=1 
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which imphes the estimate 



OO Ofc 



k=Q 1=1 





1^ — n fj — 1 



This proves the fohowing result: 

Theorem 27 Let f G (B) . Let the points {rk,e;j}^^Q belong to the interval 
[0,b] with b < 1. Then the polyharmonic cubature formula defined by 156]) with 
remainder E (/) — Jg f (x) d^ (x) — Cn (/) satisfies the following estimate 



7 Conclusions 

1. Our research may be considered as a contribution to the topic of ana- 
lytic continuation of solutions to elliptic equations (in particular, harmonic 
functions), see the discussion and references to the works of V. Avanissian, 
P. Lelong, C. Kiselman, J. Siciak, M. Jarnicki, T. du Cros, on p. 54 — 55 in 
[2], [3], [38], [15], [26], [29]. Our construction of r— analytic continuation is 
applicable to domains as annuli, strips and other domains with symmetry 
in W^. The case of the annulus has been considered in j30j, while the case 
of strip and other domains will be considered in |28| . 

2. The concept of polyharmonic Hardy spaces appears to be a new multivari- 
ate concept which differs from the existing approaches in several complex 
variables, cf. [33], [M], [IS], [H], [5], [H], [SO]. 

3. We have seen that the space of r— analytic functions on the Klein-Dirac 
quadric provides an useful setting for estimation of the remainders in In- 
terpolation and Cubature. Although the space of such functions is 1 — 1 
mapped to a Hardy space of holomorphic functions of several complex 
variables on the Lie ball, our approach has a non-trivial counterpart on 
the annulus which is not obtained from constructions, cf. [311. Our 
approach which is based on r— analytic continuation of solutions to ellip- 
tic equations (in particular, polyharmonic functions) provides non-trivial 
constructions of Hardy spaces on complexified annulus, strip and other 
symmetric domains in R'^, which are not obtained by the standard ap- 
proach to holomorphic functions in C^. 
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